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The  two  previous  technical  notes  (No.  3 and  No.  4)  described  the  use  of  the. 
data  from  tJ)C  Initial  experiment  in  development  of  an  Injury  prediction  model. 
The  present  technical  note,  whlcli  refers  to  the  discussion  ami  nutation  of  lac 
two  previous  technical  notes,  outlines  the  use  of  the  same  data  to  test  the 
reasonableness  of  proposed  head/ncck  models. 


Background  and  Assumptions 

Within  each  cell  (i,j,k,m)  of  the  initial  experiment,  there  will  be  vector 
observations 


^Ijkmn 


^ijkmnl 


^ijkmnl 


for  n “ 1,...,N  where  the  variables  yijkmnp  the  set  of  variables  describing 

dynamic  response  severity  which  were  defined  in  Technical  Note  No.  3.  For  purposes 

of  te.^tinj;  r pioposed  liead/neck  model,  it  makes  no  difference  whether  these  variables 
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are  the  variables  h or  the  transformed  variables  h. 

Por  a head/neck  model,  let  yijicmn  denote  the  vector  of  the  predicted  values 
of  yijkmn  • will  be  assumed  that  these  predicted  values  were  determined  with- 
out using  the  observed  values  in  the  initial  experiment.  Define  the  vector  of 
differences  between  observed  and  predicted  values  as 


Assume  that  yijkmn  ^ multivariate  Normal  distribution  with  mean  vector 

u. ..  and  covariance  matrix  Vi , and  that  y.  is  from  a multivariate  Normal 
_ijkmn  jkmn 

distribution  with  mean  vector  Y.  . and  covariance  matrix  V«.  If  model  predic- 

_ijknn 

tlons  were  unbiased,  then  linknn“>iikmSO  that  each  would  have  au  expected 

value  of  zero.  Thus,  a test  of  the  hypothesis  E(dijj^ny^)  “ £ will  provide  infor- 
mation about  the  usefulness  of  the  model. 

Because  yijkmn  ^ijkmn  assumed  to  be  multivariate  Normal,  it  follows 
that  the  vectors  dijj^^jj  constitute  a sample  of  size  IJKMN  from  a multivariate 
Normal  distribution.  It  should  be  noted  that  because  observations  have  been 
obtained  from  each  subject  under  all  UK  run  profiles,  there  will  be  correlation 
between  those  observations.  Fortunately,  however,  when  working  with  the  d's,  this 
poses  no  problem.  To  Illustrate  this,  consider  two  observations:  y^  for  a particular 

subject  for  a given  run  profile,  and  y2  for  the  same  subject  for  a different  run 

profile.  Let  and  corresponding  model  predictions.  Denote  the  variance 

2 

of  h^  by  and  the  correlation  between  h^  and  h2  by  p.  If  the  proposed  head/neck 
model  provides  valid  predictions,  it  is  reasonable  to  assume  that  the  variance  of 
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hj  Is  approximately  and  the  correlation  between  and  h2»  between  and  h^, 
and  between  hj^  and  ^2  is  approximately  p.  From  this  it  follows  that  the  correlation 
between  d^  » ^2  ' ^2~^2  approximately  zero,  so  that  the  vectors  dj^jj^pjj, 

may  be  regarded  as  independent.  Therefore,  if  " — » these  vectors  comprise 

a random  sample  of  IJKMN  independent  observations  from  the  same  multlvariat  3 Normal 

t 

distribution  having  zero  mean  Vector. 


Statistical  Test 

2 

In  view  of  the  above  discussion.  Hotelling's  T statistic  may  be  used  to  test 
the  hypothesis  that  E(d^jj.^)  = £.  This  statistic  Is  given  by 


- IJK>IN(d )'V“^  (d ) 

where  d is  the  (P  x 1)  mean  vector 
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and  V is  the  (P  x P)  sample  covariance  matrix 


{ EE  ■‘Utanp»lik»nq  - 


IJKMN  ^ ijkmn 
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The  hypothesis  may  be  tested  at  an  a significance  level  by  comparing  the  observed 
2 

T value  to  the  value 


T.  ‘ - IJKMN-1  F„  „ . 

“ TjkwT:? 

where  Fp  q represents  the  upper  a point  of  an  F distribution  with  P and 

IJKMN-P  degrees  of  freedom. 
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If  T >T(j  , the  hypothesis  that  ••  0 Is  rejected,  Indicating  that 

the  proposed  model  is  inadequate.  In  this  situation,  univariate  tests  on  each 

^IJkmnp  pl-oCs  of  these  variables  may  reveal  the  source  of  the  model  Inadequacy. 

2 2 

If,  on  the  other  hand,  T £T  this  docs  not  necessarily  mean  that  the  model  is 

^ 4 

completely  adequate,  but  only  that  there  is  not  enough  evidence  in  the  data  to 
disprove  it.  If  a model  does  pass  this  test,  the  variability  of  the  predicted 
values  should  be  examined  relative  to  the  inherent  variability  in  the  observed 
dynamic  re'*ponse  variables. 
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